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Summary 


y This  report  reviews  initial  progress  on  a new  program  of  investigation 
of  scintillation  and  coherence  effects  for  a laser-illuminated  , noncoopera- 
tive target  as  viewed  through  atmospheric  turbulence.  The  application  of 
the  effort  is  in  the  prediction  of  turbulence  effects  on  the  operation  of 
coherent  optical  adaptive  transmitter  (COAT)  systems. 

Significant  analytical  progress  is  reviewed,  including  results  for  the 
mutual  coherence  function,  variance  and  covariance  of  irradiance,  spectra, 
and  statistics  for  a coherently  illuminated  diffuse  target.  The  approach 
utilizes  the  extended  Huygens-Fresnel  principle,  and  includes  turbulence 
effects  in  both  the  path  from  the  transmitter  to  the  target  and  that  back 
to  the  receiver.  It  is  found  that  there  are  three  pertinent  covariance 
sc  * i ~s  and  six  possible  parameter  realms,  and  that  the  normalized  variance 
will  be  unity  except  in  those  cases  when  the  target  spot  is  sufficiently 
c5nall  as  to  constitute  a quasi-point-source.  This  parameter-realm  view- 
point is  further  explored  in  relation  to  a variety  of  possible  sources 
operating  through  turbulence.  The  analysis  is  partially  extended  to  a 
more  complex  target  and  to  the  incoherent  case,  and  future  analytical  tasks 
and  applications  to  real  adaptive  systems  are  outlined. 

The  establishment  of  an  experimental  field  facility  is  also  described, 
which  will  be  capable  of  measuring  all  pertinent  quantities  at  both  visible 
and  middle-infrared  wavelengths./  Preliminary  experimental  results  are 
presented. 

The  results  of  this  work  should  aid  materially  in  the  understanding  of 
the  dynamic  behavior  of  adaptive  laser  sources  operating  in  the  presence  of 
significant  atmospheric  turbulence. 
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I.  Introduction 

This  report  reviews  a new  and  ongoing  investigation  of  scintillation 
effects  arising  from  the  illumination  of  noncooperative  targets  with  laser 
radiation  in  the  presence  of  atmospheric  turbulence.  The  application  of 
interest  is  an  understanding  of  the  performance  of  adaptive  optical  transmit- 
ter systems  in  the  presence  of  target  and  turbulence-induced  speckle  and 
other  scintillation  structure. 

During  the  initial  period  of  research,  a theoretical  understanding  of 
turbulence  effects  was  achieved  for  diffuse  targets  with  coherent  and 
incoherent  (highly-multimode)  illumination.  Progress  was  also  made  on 
the  problem  o^  target  structure  (glints)  ; and  on  scintillations  from  gen- 
eralized sources,  including  parameter  realms  and  covariance  scales.  These 
theoretical  results  are  discussed  in  detail  in  Section  II. 

During  the  same  period,  an  experimental  program  was  developed  for 
visible  and  infrared  (3. 5-3.8  micron)  wavelengths.  Preliminary  field 
experiments  were  conducted  at  6328$  and  a well-instrumented  field  facility 
established  at  4880&.  Concurrently,  a pulsed  infrared  laser  and  receiver 
system  was  designed  and  fabrication  initiated.  The  experimental  efforts 
are  discussed  ::n  Section  III. 

II.  Theoretical  Description 

Significant  progress  has  been  made  in  the  physical  and  analytical 
understanding  of  turbulence  effects  cn  the  dynamics  of  target-reflected 
radiation.  The  important  quantities,  relating  ultimately  to  the  performance 
of  an  adaptive  transmitter  system,  are  the  variance  of  irradiance  (o^) , 
covariance  of  irradiance  [C  (p)],  mutual  coherence  function  [T(p)], 
probability  distribution  function  of  irradiance,  and  spatial  and  temporal 
power  spectrum.  These  quantities  are  of  interest  in  the  plane  of  the  active 
laser  transceiver,  and  include  turbulence  effects  on  both  the  illuminating 
radiation  from  the  transmitter  to  the  target  and  on  the  scattered  radiation 
over  the  return  path  (Fig.  1). 

The  treatment  to  be  given  below  utilizes  primarily  the  extended  Huygens- 
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Figure  1.  Illuminator,  Target,  and  Receiver  Configuration 
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Fresnel  principle  as  applied  to  a turbulent  path.  * Until  otherwise 
stated,  the  laser  is  assumed  to  be  a coherent  (TEMAA)  source,  collimated 
or  focused,  with  a perfectly  diffuse  target.  In  all  cases,  we  attempt  to 
show  clearly  the  assumptions  and  approximations  that  are  made  and  to  dis- 
cuss their  implications.  An  important  feature  is  the  definition  of  dis- 
tinct (asymptotic)  parameter  conditions  applying  to  any  given  configuration; 
there  are  generally  six  such  conditions,  representing  the  possible  permuta- 
tions of  inequalities  between  the  three  pertinent  parameters:  Fresnel  zone 

% 

size  (L/k)  , coherence  radius  (o  ) and  speckle  size  in  the  absence  of 

o 

turbulence.  Each  such  condition  will  in  general  carry  a distinct  physical 
and  analytical  interpretation.  It  will  be  seen  that  cases  of  strong  scin- 
tillations ("saturat ion**  or  multiple  scattering)  are  included  In  these  con- 
ditions, so  that  the  treatment  is  general. 

We  first  consider  (Sec.  II-A)  the  most  common  situation,  i.e.,  that  in 
which  the  primary  effect  of  the  turbulence  on  the  reflected  radiation 
arises  through  the  perturbation  of  the  phase  term  in  the  associated  Green* s 
function.  The  implications  of  this  assumption  in  terms  of  the  field  sta- 
tistics are  also  explained.  A review  of  some  of  the  developments  of  this 

3 

section  appeared  in  the  preceding  technical  report  but  without  a full 
interpretation  of  the  assumptions. 

We  then  generalize  the  development  (Sec.  II-B)  to  include  the  effects 
of  the  amplitude  perturbation  term  in  the  Green* s function.  The  implica- 
tions are  again  explored  in  detail.  In  Sec.  II-C,  we  treat  the  first-order 
statistics  of  irradiance  for  a target  containing  one  or  more  glints. 
Incoherent  illumination  is  then  discussed  in  Sec.  II-D. 

In  Sec  II-E,  we  relate  the  current  development  to  an  alternative, 
earlier  treatment  of  a related  problem  and  show  that  the  two  approaches 
are  basically  complementary.  We  also  consider  the  possible  parameter  realms 
for  a variety  of  source  types,  with  attendant  predictions  of  strengths  and 
spatial  scales  of  scintillations. 

Finally,  in  II-F,  we  outline  further  topics  and  extensions  of  interest 
and  related  efforts  to  be  undertaken. 

1.  R.F.Lutomirski  and  H.T. Yura, "Propagation  of  a Finite  Optical  Beam  in  an 
Inhomogeneous  Medium",  Applied  Optics,J_0,  1652,  July  1971. 

2.  H. Yura, "Mutual  Coherence  Function  of  a Finite  Cross  Section  Optical  Beam 
Propagating  in  a Turbulent  Medium", Applied  Optics,  H_,  1399,  June  1972. 

3.  J.R.Kerr,  et  al,  "Propagation  of  Multiwavelength  Laser  Radiation  through 
Atmospheric  Turbulence",  RADC  Technical  Report,  August  1975. 
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II. A, Basic  Irradiance  Statistics  and  Mutual  Coherence  Function 

Previous  work  on  speckle  statistics  has  primarily  been  concentrated 

on  the  nature  and  statistics  of  the  target  surface,  propagation  of  the 

speckle  field  without  turbulence,  and  effects  of  speckle  on  image  quality. 

Speckle  propagation  through  turbulence  has  been  considered  over  a vertical 

4 

path  for  the  purposes  of  speckle  interferometry. 

In  the  present  section  an  analysis  is  given  of  the  first  and  second 
order  statistics  of  the  received  intensity  (irradiance)  after  scattering 
from  a diffuse  target.  The  treatment  is  based  on  the  extended  Huygens- 
Fresnel  formulation  and  includes  the  effects  of  the  turbulent  atmosphere 
on  the  laser  beam  as  it  propagates  to  the  target  and  on  the  speckle  as  it 
propagates  back  to  the  receiver,  Formulations  arc  given  for  both  the  focused 
and  collimated  cases.  The  analysis  also  includes  the  mutual  coherence  func^ 


tion  (MCF) , 

The  source,  target,  and  receiver  configuration  is  shown  in  Figure  1. 
The  present  analysis  is  confined  to  the  case  of  a TEMqq  laser  illuminator. 
The  source  and  target  are  assumed  to  be  much  smaller  than  the  path  length 
(L) , and  the  distance  between  the  receiver  and  source  is  greater  than  the 
source  size  and  much  smaller  than  the  path  length.  These  geometric  con- 
ditions confine  the  problem  to  small  angles  and  ensure  that  the  outgoing 
and  returning  radiation  experience  independent  turbulence  regions;  the 
latter  limitation  is  thought  to  be  Inessential  owing  to  the  diffuse  target 
characteristics. 

1.  Mean  Irradiance  at  Receiver 

To  find  the  mean  irradiance,  we  need  no  assumptions  other  than 
that  of  a diffuse  target. 

We  write  the  source  amplitude  distribution  as 


U (r)  = U exp 
o o 


2T 


(1) 


and  F are  the  characteristic  beam  radius  and  focal  length  respec- 

A.  M.  Elbaum,  et  al,  "Laser  Correlography«  Transmission  of  High-Resolution 
Object  Signatures  through  the  Turbulent  Atmosphere",  Riverside  Research 
Institute,  Technical  Report  T-l/306-3-11 , October  31,  197A. 
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tively.  The  field  at  the  target  is  written  from  the  extended  Huygens- 

1 2 

Fresnel  principle  * as 


U(p; 


k ikL  f 

2? Il"  J Uo(r)  exp 


ik  1 p - r 1 
2L 


+ <l»i  (P,r) 


dr 


(2) 


where  describes  the  effects  of  the  random  medium  on  the  propagation 
of  a spherical  wave.  Combining  Eqs.  (1)  and  (2),  we  have 


o 


+ V ) (P,r) 


dr 


(3) 


In  particular,  this  applies  to  the  special  cases  of  a focused  (L  = F) 
or  collimated  (F  -+■  ®)  beam  respectively. 

The  field  at  the  receiver  is  written  by  reapplying  the  Huygens- 
Fresnel  principle  to  the  field  at  the  target: 


ik  [L  + P2; 

] 

..  -v  ke  L 21. 

u<P)  - - 2VlL  - 

-l  U'(p)  exp  (p2  - 2p-p)  + dP  ^ 

where  U'(7)  is  the  field  solution  after  reflection  from  the  target,  and 
\p2  represents  the  turbulence  effect  from  the  target  to  the  receiver.  The 
mean  intensity  at  the  Receiver  is  then 

<l(p)>  = <|U(p)|2>  = Jj  dpjdpj  <U' (pi)U'*(p2)> 

* exP  [§  ((pf  - Pi)  - 2p.  (P!-P2))]  < 

• < exp  [*2(p,Fi)  + t^*(p,P2)]  > 
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Through  the  assumption  of  a diffuse  target,  the  reflected  beam 
suffers  a random  phase  delay  from  point-to-point  over  the  target,  so  that 


<u' (Pi)U’*(p2)>  = <I(Pl)>  6 (p j — p 2 ) 


(6) 


Using  this  in  Eq.  (5),  the  mean  intensity  becomes 

2 C 

<I(p)>  = (JtTl)  J dpl  <lL'(Pl)l'>  < exp  j^2(p,Pi)  + V<2*(p,Pi)J> 

where  the  mean  exponential  term  is  unity  from  considerations  of  energy 
2 

conservation.  The  resultant  mean  intensity  at  the  receiver  is  then 
simply 


(7) 


<Kp)> 


(■«7 


dp  <|U(p)|2> 


(8) 


Tp  complete  the  solution,  we  use  Eq,  (3)  with  Eq.  (8).  We  note 
that  the  structure  function  gives  us  (r  « | r j -r"2  I ) 


<exp  ^i(p,ri)  + ih*(p,r2)J 

For  the  focused  beam,  we  then  have 


/ r\ 

'Vo) 


5/3 


> = 


(9) 


'u(?)l2>  ■ (sl)  2 -//  <*ndr2 

/ \ 5/3~ 

-fc)  J 


exp 


•f  + 


2a  2 
o 


ik  — — — . 

"l"  P ’ (r l-r2> 


(10) 


Carrying  out  the  integration  indicated  in  Eq.  (8),  involving  the  Fourier- 
Bessel  integral,  we  have  finally 

)20’2  ? 
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The  result  for  the  collimated  beam  is  identical,  and  in  fact  could 
b«  deduced  for  an  arbitrary  beam  focus  (Eq.  (1))  through  conservation  of 
energy: 

<I(p)>  = (ffe)2/ d7  <m\*> 


0 


Thus  the  mean  intensity  at  the  receiver  (illuminator)  plane  is  uniform  and 
independent  of  turbulence  level. 

2.  Correlation  Function  of  Irradiance 

In  order  to  calculate  the  correlation  function  or  covariance 
of  irradiance,  we  assume  for  the  present  section  that  the  perturbation  Green’s 
function  (wave  structure  function)  is  dominated  by  the  phase  perturbation 
(phase  structure  function),’*  This  will  be  true  for  many  cases  of  interest 
and  will  be  relaxed  in  subsequent  sections,  where  the  actual  implications 
of  the  assumption  are  pointed  out. 

The  correlation  function  of  the  intensity  at  receiver  points  pj 
and  P2  is  given  by 


BI<'Pl»P2)  * 


*1  (Pl)I2^P2^>  " <u(Pl)'J*(Pi)U(p2)U*(p2)>  (13) 


Utilizing  the  extended  Huygens-Fresnel  principle,  the  correlation  function 
can  be  expressed  as 


Bj(Pl »P2> 


dp jdp2dp 3dp4  <U(p1)U*(p2)U(p3)U*(p4)> 


exp[  ik(Rj  R2  + r3  _ R^ ) ] H(pj  »P2,P.3,Pil;p1  ,p2)  (14) 


5,  R.  S.  Lawrence  and  J,  W,  Strohbehn,  "A  Survey  of  Clean-Air  Propagation 
Effects  Relevant  to  Optical  Communications",  Proc,  IEEE,  58,  1523, 
October  1970, 
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where  H is  the  fourth  order  mutual  coherence  function  given  by 


H - < exp[iHP],Pi)  + <P*(p2,Pi)  + ^(p3,p2)  + ^*(p4,P2)i  > 


(15) 


and 


Rl  = 
R2  = 

*3  = 
RU  = 


IP1 
I Pi 
IP2 

IP2 


Pi 

P 2 
1 — . 
P3  I 

P4  I 


Under  the  assumption  of  dominant  phase  perturbations, 

H = < exp[i<  (P]  ,Pi)  - i^> (o2 ,Pi ) + i<HP 3 » p2 ) - > (16) 

After  reflection  from  the  diffuse  target,  the  fields  are  Gaussian 
and  spatially  incoherent.  Therefore,  the  fields  at  the  target  can  be 
expressed  as 


<U  (p  j )U"  (p2)U  (p  3 )U"  (o,. ) > = <U(p1)U*(p2)>  <U (p*3 ) U*  (o^ ) > 
+ <U(P])U*(p4)>  <U*(p2)U(p3)> 

= <I(Pl)>  <I(p3)>  «(pj-p2)  0(^3-^) 


+ <I(Pi)>  <I(P3)>  6(Pl-p4)  6 (p 3~p2 ) 

Utilizing  (17)  and  (14)  the  correlation  function  can  be  expressed  as 


(17) 


BI(Pl»P2)  = (y^)  jjjT  dP2  do4  <1(p2)>  <I(P4)> 


fea  j] 


T _ _ 


dp2  dp4  <I(p2)>  <I(p4)>  exp  ~^-(p2_p4)  (pl 


H(p2,p2»P4>P4>Pl>P2) 


(18) 
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where  use  has  been  made  of 


R>-R2  = 2L 


R3-R-  ~ 2L 


P ] P 2 2 (p ] “p2 ) * p ] 


P32  P4 ^ - 2 (p 3 - p ^ ) • p 


The  fourth  order  mutual  coherence  function  (see  Appendix  A)  in  Eq.  (18) 

, v 2,6,7 

is  given  by 


*’(P2  >P2  >P4  »p4  *P]  > P2  ) “ H (Pj  , p2  ,P  3 ,P4  ;pj  ,p2  , P3  ,P4) 

Pi  = P2 

P3  = P4 
P3  = Pi 
P 4 = P2 


|(D12  “ D13  + D14  + D23  “ °24  + D34) 

= e 

P]  = P2 
P3  = P4 

P3  = Pi 

P4  = P2 

where  the  wave  phase  structure  functions  D are  given  by 


(19) 


6..  D.  L.  Fried,  "Effects  of  Atmospheric  Turbulence  on  Static  and  Tracking 
Optical  Heterodyne  Receivers",  Optical  Science  Consultants,  Technical 
Report  TR-027,  August  1971. 

7.  D.  L.  Fried,  "Atmospheric  Modulation  Noise  in  an  Optical  Heterodyne 
Receiver",  IEEE  Trans,  on  Quantum  Electronics,  QE-3,  213,  June  1967. 
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I 
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rorwwwwpf 


= p 5/3 


I Ip^-Pi 1 8/3-  Ip.-pJ873! 

A. 


(Pj-Pi)  ~ (P.-P.)| 


(20) 


where 


P„  = (0.545  C 2Lk2) 
o n 


-3/5 


scale  and 


is  the  turbul jnce-induced  coherence 


2/3 

cnz  = Structure  constant  of  index  of  refraction  (nf  ). 
Using  this  in  (18)  and  making  the  change  of  variables 


and 


P - P2  - P4  , p = pi  - p2 


2R  = p2  + P4 


and  recognizing  that  the  first  term  in  Q8)  equals  <I(p1)><I(p2)>,  the 
covariance  for  the  focused  case  is  given  by 

5/3 

(i  / O » O / T\  \ 

Cl(?,,?2)  - (,--L)  (L-)  |Uj 


'a  2\  2 -2^> 

_o__ 

. 2 ~j 


e V3C 


rR-ty*L 


ik  — — 

P • p 


7/ 


r ] r2^r 1 ^r2  exP 


(r . 2 » 2\  / ^/3  5/3 

(r 1 + r2  ) (rj  + r2  ) 


4a  2 
0 


5/3 


5/3 


p5/3  _ I IP8/3  ~ P6/3! 

2 Ip  - p| 


_ I lp8/3  - 08/3| 

Ip  + p|  J 


-11- 


(21) 


* 


The  covariance  for  the  collimated  case  is  obtained  from  (21)  by  replacing 
the  dr^dr2  integration  by 


r2dr 


jdr2 


(r-2 


1 + r20 


|p8/3  - pfl/3|  2 Ip8/3  - p 8 / 3 ! 

•y 

I p - p I I p + p I 


(22) 


In  order  to  further  reduce  the  number  of  integrations,  the  zero 
order  Bessel  functions  must  be  expanded  to  functionally  separate  the  R and 
p dependence.  This  can  be  accomplished  by  utilizing  the  following  iden- 
tities : 


= X 

in=0 


c 

m 


m 

(+1) 


ID 


(l  rJR)  Jm  (l  rh)  COS  " 4 


where 


and 


e=2 

ra 


for  m ^ 0. 


The  integration  is  then  given  by 


(23) 
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J Jo  (t  r‘|R  + l1)^  (l  r2'R  ' I1)  d0R 
0 

■ 2>  rn  <-1)" '» (r  r>R)  Jm  (t  r^R) 

nr  *U 

• (l-  r»!)  \ (c  r4) 


Using  the  Fourier-Bessel  integral, 


(24) 


Jrjdnffr,)  / RdR  ^ r 
o o 


and  the  covariance  for  the  focused  case  becomes 


= f (r2)for  r2^0 

2 

_2/L\  f(o)  for  r2=0. 

(25) 


focused 


• exp 
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V 1 

rz  j 
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ik  — — 2 

- p • 0 

p 5/  3 _ I |p8/3  - P8/3|  1 J p ® / 3 _ p8/b7| 

p 5/3 

— 0 

lp  - p|  2 |p  + p| 

where  lse  has  been  made  of 


(26 ) 


Z 

m=0 


m 

(-1)  e 


m 


J 2 
m 


(x)  = Jq  (2x) 


(27) 
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The  variance  is  then  given  by 


and  6 is  the  Kronecker  delta. 
r 2 


and  the  normalized  variance 


is  unity  (Eq.  12) . 

The  same  technique  yields  the  covariance  for  the  collimated  case 
It  is  given  by 


r2^r2 


i2r 
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>/: 
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/dir  j 

J o 


L ra 


p) 


6 

r2 


] 


p 5/ 3 _ I [p8/3  - p 8 / 3 | 1 | p 8 / 3 _ p8/3j 

Ip  - p\  2 Ip  + p] 


and  the  normalized  variance  is  again  unity. 


(29) 
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The  above  results  are  based  on  the  assumptions  of  a diffuse 

target  and  phase  perturbations  dominating  the  turbulence  effects.  The 

resulting  threefold  integrations  will  be  numerically  evaluated  in  future 

work.  Physical  interpretation  is  not  difficult  but  will  be  clarified 

in  the  discussion  of  further  simplifications  below.  It  suffices  to  point 

out  here  that  the  covariance  fundamentally  involves  two  scales:  the 

turbulence-induced  coherence  scale  (PQ)  and  the  speckle  scale  in  the 

absence  of  turbulence  (a  and  L/ka  for  the  focused  and  collimated  case 

o o 

respectively) . 

Also,  we  note  that  within  the  present  assumptions,  the  normalized 
variance  of  intensity  is  always  unity  independent  of  turbulence  strength. 
This  tgrees  with  a physical  model  of  identical-frequenc}*,  randomly 
phasea  oscillators  summed  to  Tepresent  any  given  point  in  the  receiver 
field:  the  model  applies  regardless  of  whether  target  speckle  (aQ  or 

L/ka^)  or  Matmospheric  speckle*'  (pQ)  dominates. 

3.  Mutual  Coherence  Function 

The  mutual  coherence  function  (MCF)  may  be  very  important  in 
analyzing  the  operation  of  a coherent  optical  adaptive  system,  and  can 
be  readily  derived  given  the  assumption  of  a diffuse  target  but  without 
assuming  dominance  of  the  phase  perturbation  term.  We  write 


r(pi ,P2) 


d p j d p 2 


< U(Pl)U*(P2) 


> 


exp 


(P i .P\) 


exp 


JjMPl.Pl) 


+ ^2*  (P 2 


» P2 )] 


(30) 


where  R](p]»Pi),  ^2(P2>P2)  are 'the  distances  from  pj  to  pj  and  p2 
to  p2  respectively. 


By  the  Fresnel  approximation 

Rl(p],Pl)  - R2(P2,P2)  = ELl±1±?LA  - • Pl  ~ P2  • P~2 

2L  L 


(31) 
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Finally,  from  (30)  and  (31), 


Since  the  wave  is  incoherent  after  reflection  from  the  diffuse 
target,  the  coherence  function  at  that  plane  can  again  be  represented  by 
the  Dirac  delta  function  as  given  in  (6).  Using  this  in  (32),  r(pltp2)  can 
be  simplified  to 


lk  Cp i -p| ) 


2L 


I 

/ 


dp  <I(p)> 


cxp  ik  (P1-p2).p  | 


(33) 


In  the  absence  of  turbulence,  this  equation  is  entirely  identical  to  the 
Van  Cittert-Zernike  theorem  of  coherence  theory,9  which  is  identical  to  a 

result  obtained  by  Goodman  for  the  mutual  coherence  function  of  a pulsed 
optical  radar.10 

To  complete  the  solution,  we  utilize  the  mean  intensity  at  the 
target.  For  the  focused  case: 


<!(?»- (i)2  l».l 2 “4- 


We  thus  have 


9.  M.  Born  and  E.  Wolf,  Principles  of  Optics.  Perg&mon  Press,  New  York. 
1975. 

10.  J.  W.  Goodman,  "Some  Effects  of  Target-Induced  Scintillation  on 
Optical  Radar  Performance",  Proc.  IEEE,  53,  1688,  1965. 
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From  the  Fourier-Bessel  integral  formula. 
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Wo  Vl  rp)  Jo(lPC)  d3 
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Equation  (34)  can  then  be  simplified  and  it  becomes 


r<P”P2)  ‘ £(t)  |Uo|2  V e 


J£i_.  2/£ 

4a  2 


5/3 


ik 

2L 


(p?”P^) 


focused 


_ _ 2 P- 
4a  2 Vp' 


5/3 


ik 

2L 


(Pif-Pl) 


= <I(p)>  e 


(36) 


Using  <I(p)>  for  the  collimated  case  in  (33)  and  simplifying  yields 


r Cp  1 * P2 ) 55  e 

Collimated 

These  results  for  the  MCF  will  be  used  further  below. 

It  inay  be  noted  that  the  MCF's  of  Eq.  (36)  and  (37)  imply  a "white" 
or  constant  spatial  power  spectrum  for  the  (complex)  amplitude.  This  is 
of  course  an  idealization  resulting  from  assuming  delta-function  rather 
than  wavelength-sized  phase  correlation  for  the  field  upon  reflection  off 
the  target.  The  more  interesting  spectrum,  however,  is  that  of  the  irra- 
diance,  as  discussed  below, 

4.  Probability  Distribution 

In  order  to  formulate  the  probability  distribution  for  the  scin- 
tillating energy  at  the  receiver,  we  evaluate  the  nth  moment  of  the  intensity 
in  terms  of  <I(p)>.  We  again  assume  that  the  phase  perturbations  are  the 
dominant  turbulence  effect,  and  write  the  second  moment  as 

• <U(p,)u*(F2)U(?3)UO(^)>  exp  ("‘-is 

+ pi  - p?  - 2^,-^) . p-2(p3.p4)  .-yi 


• <«P  ji(*(Pi>p)  - ♦(»„,)  + *;p3,p)  . +(-  -;j 

Since  the  fieldd  after  reflection  from  the  diffuse  target  are  jointly 
Gaussian  and  uncorrelated, 


(39) 
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Using  (40)  in  (38)  yields 


4 

<l2(P):  = 2 ('2wlJ  ff  dp2dp4  <I(P2)>  <I(P4)> 

Since 

<I(p)>  = (ife)  / dp2  <i(p2)> 

the  second  moment  becomes 

<1  (p)>  = 2<I(p)>2 


(41) 


(42) 


Similarly,  it  may  be  shown  that  the  nC^  moment  is  given  by 

<ln(p)>  = N!  <I(p)>N  (43) 
The  probability  density  function  for  the  intensity  therefore  is  exponential 
and 


J>j(a) 


_e_ 

<I> 


<I> 


(44) 


where  a is  greater  than  or  equal  to  zero.  It  is  thus  concluded  that  the 
field  and  amplitude  at  the  receiver  are  normally  and  Rayleigh  distributed 
respectively,  given  the  assumption  that  phase  perturbations  dominate  the 
turbulence  effects. 


5.  Simplification  for  Weak  and  Strong  Turbulence 

Since,  as  shown  above,  the  fields  at  the  receiver  are  gaussian 
and  spatially  "white",  it  is  tempting  to  assume  that  the  receiver  fields 
are  also  jointly  gaussian.  This  turns  out  to  be  a good  approximation  in 
many  situations,  and  in  this  section  the  implications  and  conditions  for 
validity  of  this  added  assumption  are  explored.  This  leads  to  a simple, 
straightforward  interpretation  of  the  terms  in  the  covariance  of  intensity. 

The  jointly  Gaussian  assumption  yields 

= <U(p1)U':(pi)><U(p2)bA(p2)>  + <lUPl)U*(p2)><UA(pOU(p2)> 

_ _ 2 

= <l(p])><I (P2 )>  + I ^ (p i » p2 ) I (45) 

It  follows  that  the  covariance  of  intensity  is  given  by 

CI(Pl»P2>  = bj  (p  1 *P2)  - <1  (P])><'I  (P2>>  = ! f (p  ] > P2 ) I (46) 


Finally,  utilizing  the  mutual  coherence  result  (Eq.  36)  the 
normalized  covariance  function  of  irradiance  for  the  focused  case  can  thus 
be  written: 


focused 


where  the  normalized  variance  is  unity  as  before. 

For  the  collimated  case,  the  same  variance  is  obtained,  and  the  normalized 
covariance  is 


-A 

CT  (p)  = e 
N 

collimated 


(48) 
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The  covariance  scale  lengths  are  obvious  from  these  results. 

Either  the  "atmospheric  speckle"  (p^)  or  the  "target  speckle"  (speckle 
in  the  absence  of  turbulence)  will  predominate,  depending  upon  which  is 
smaller  (strong  and  weak  turbulence  respectively).  We  point  out  in  passing 
that  a third  covariance  scale  (L/k)^  may  also  enter,  but  this  scale  is  lost 
within  the  present  assumption  of  dominant  phase  perturbations  (see  Sec- 
tions  II-B  and  II-E  below) , 

We  note  that  the  spatial  power  spectrum  of  inadiance  may  be 
readily  obtained  by  transforming  Eqs, (47,48),  However,  a more  important 
quantity  in  the  operation  of  e.g.  an  adaptive  optical  system  may  be  the 
temporal  spectrum.  This  spectrum,  which  will  be  derived  in  a later  sec- 
tion, depends  only  on  the  atmospheric  speckle  term;  the  target  speckle 
field  will  not  translate  with  the  transverse  wind. 

We  now  explore  the  conditions  for  validity  of  the  jointly  gaussian 
assumption.  The  simple  multiplicative  terms  for  the  covariance  scales  in 
Eqs.  (47,  48)  are  replaced  bv  a more  complicated  interrelationship  in 
Eqs.  (26,29),  It  may  therefore  be  surmised  that  the  jointly  gaussian 
assumption  is  valid  under  conditions  of  weak  and  strong  turbulence,  when 
target  and  atmospheric  speckle  terms  respectively  predominate,  but  that 
the  jointly  gaussian  description  is  not  correct  in  the  range  of  inter- 
mediate turbulence  effects  when  both  scales  are  important  and  interact, 

We  now  show  that  this  is  indeed  the  case. 

Weak  Turbulence 

For  the  weak  turbulence  case,  p »/hJk  and  the  term 


5/3 

P 


8/3  8/3, 

P - "P  I 

I p~p 


1 

2 


(49) 

in  (26)  and  (29),  The  covariances  then  become  identical  to  those  derived 
using  the  jointly  gaussian  assumption. 


Strong  Turbulence 

In  the  strong  turbulence  case,  po«/L/k,  which  corresponds  to 
multiple  scattering,  or  "saturation  of  scintillations"  for  a point  source. 
Let  us  consider  Eq.  (26)  with  p^-K),  The  only  interesting  range  of  the 
argument  (p)  is  0 <p<p  . The  Bessel  term  O^pJ  in  polar  coordinates)  is 
appreciable  only  for  p^O,  i,e,  P>>Po  or  p»p;  and  because  of  the  latter 
condition  the  final  bracket  in  the  equation  is  zero.  Hence  the  condition 
(49)  is  again  obtained,  and  the  covariances  again  become  identical  to 
those  derived  using  the  jointly  gaussian  assumption.  The  atmosphere  has 
"decoherentized"  the  field  in  a manner  similar  to  that  of  the  diffuse 
reflector  itself, 

6.  Time-Lagged  Covariance  and  Temporal  Spectrum  of  Scintillations 

For  this  development  we  assume  that  the  fields  are  Jointly  gaussian 
at  the  receiver  and  consequently 

Vpi.Fj.t'  - |r(p1J2>T)|2  (50) 

Using  the  extended  Huygens-Fresnel  principle 

ik(pi2-p22) 
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- 2p2 

p2)J  < exp  Hji 

2(pl»Pl»°)  + ^2*(P2»P2>T) 

Due  to  the  diffuse  target 

<U  (P1,0)U*  v'p2,t)>  = (Pi,0)U*(fT2,T)>6(p1-p2) 


(52) 


Using  (52)  in  (51)  and  utilizing  the  extended  Huygens-Fresnel  principle 
to  express  the  fields  incident  on  the  target,  the  time  delayed  mutual 
coherence  function  at  the  receiver  becomes  for  the  focused  case 
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(53) 


The  first  ensemble  average  in  (53)  corresponds  to  the  time  delayed 
mutual  coherence  function  for  spherical  waves  originating  at  two  points 
r^  and  r 2 in  the  transmitter  plane  and  propagating  to  a single  point  pj 
in  the  target  plane.  The  second  ensemble  average  corresponds  to  the 
time  delayed  mutual  coherence  function  for  a spherical  wave  originating 
at  the  point  pj  in  the  target  plane  and  propagating  to  two  points  p^  and 
P2  in  the  receiver  plane.  Performing  the  integrations  in  (53)  it  becomes 


f(P.0,t)F(0,p,t)  (54) 


where  the  time  delayed  mutual  coherence  functions  F(p,0,t)  and  F(0,p,t) 
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can  be  obtained  from  the  mutual  coherence  function  for  a spherical  wave 
by  invoking  Taylor's  Hypothesis.5 

This  mutual  coherence  function  is  given  by2 


I:(r»p)  = exp  -2.91  Lk2 


f V ^ lwp  + (l-w)r  | dw 


(55) 


where  r is  the  transmitter  aperture  vector  r1-r2,p  is  the  target  aperture 
vector  Pi-p2  and  w is  the  distance  from  the  source  to  the  field  point 
normalized  by  the  total  path  length  L.  For  the  uniform  turbulence  case, 
we  note  that  k2LC^2'\/po  5/?.  The  time  delayed  mutual  coherence  function11 
can  be  obtained  from  (55)  by  replacing  wp"  by  wp  - V(w)t,  where  V(w)  is 
the  transverse  wind : 


p(r»P»T)  = exp 


-2.91  Lk2 


J cn2(w)|wp  - V (w) t + (l-w)r|5/3  dw 


and  consequently 


(56) 
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exp 
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o J 


and 


f(0,p,t)  = exp 
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(58) 


where  z = (1-w) . 

The  normalized,  time  delayed  covariance  function  for  the  focused 
case  and  an  arbitrary  distribution  of  Cn2  along  the  path  is  thus  given  by 
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11.  R.S.  Lawrence,  G.R.Ocho,  and  S.  F.  Clifford,  "Use  of  Scintillations  to 
Measure  Average  Wind  Across  a Light  Bean",  Applied  Optics, 11,  239 
February  1972.  — * 
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For  the  collimated  case,  a further  multiplicative  factor 


completes  the  expression. 


The  temporal  socctral  density  of  irradiance  can  be  obtained  from 
(59).  Letting  p = 0 and  assuming  uniform  turbulence  and  crosswind,  (59) 
becomes 
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Taking  the  Fourier  transform  of  (60)  yields  the  spectral  density 
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where 


|v| 


The  normalized  spectral  density  is  plotted  /ersus  the  parameter  (x)  in 
Figure  2. 
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II. B. The  General  Covariance  Function 


In  the  previous  section,  the  mean  irradiance  and  mutual  coherence 
function  were  derived  with  the  simple  assumption  of  a diffuse  target.  How- 
ever, the  analysis  of  the  intensity  covariance  was  based  on  the  added 
assumption  of  neglecting  the  amplitude  term  in  the  perturbation  Green's 
function.  This  latter  assumption  is  not  always  valid,  as  can  be  seen  by 
considering  the  case  where  the  illuminated  target  spot  in  small:  the  scin- 
tillations at  the  receiver  should  then  approximate  those  of  point  source, 
with  log  normal  statistics  and  a covariance  scale  on  the  order  (L/k)*5. 

In  this  section  we  include  the  amplitude  perturbation  term  and  examine 
the  effect  on  the  probability  distribution  of  the  irradiance.  We  then 
derive  the  general  covariance  function,  with  physical  interpretation.  The 
actual  conditions  for  validity  of  the  simpler  result  of  Sec.  II. A will  also 
be  discussed,  with  further  elaboration  in  Sec.  II. E. 

1.  Moments  and  Probability  Distribution  of  Irradiance 

The  second  moment  of  irradiance  can  be  written  from  Eqs.  (14)  and 
(17)  and  the  mutual  coherence  function  (15).  We  now  generalize  the  latter 
i-o  include  amplitude  perturbation  terms  (Appendix  A),  resulting  in 

<I2(P)>  = 2(2^)  JJ  dp2dp4  <I(P2)><I(P4)>  e4Cx  ^ I P 2 ~ I ) (61) 

where  is  the  point-source  log  amplitude  covariance  function  given  in  the 
first-order  theory5* 12 * 13  by 

» L 

yfo-ej)  - 4»vJ  f u»(u)Jo(l^£i.lus)  Sln2[if|gzsl]dsdu  (62) 

0 0 

and  $ is  the  Kolmogorov  spectrum5 


$(u)  = 0.033  C 2u"n/3 
n 


(63) 


the  present  that  V*  whlch  w111  be  t™e5  for  both  weak 

12’  M on  TMtMSkli’  ^e-Effects  °f  the  Turbulent  Atmosphere  on  Wave  Pronav; 
I,  =r^j*  National  Technical  Information  Service  (0TT-68-5O464) , 1971. P * 

* ee  an<*  Harp,  "Weak  Scattering  in  Random  Media  with  Add!  i 

cations  to  Remote  Probing",  Proc.  IEEE,  57,  375,  April  1969.  PP 
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and  saturated  scintillations  (weak  and  strong  turbulence  scattering  respect- 
ively) . It  will  not  be  true  for  the  intermediate  case,  however,  and  we 
will  lift  the  restriction  in  calculating  the  covariance  below.  Also,  we 

point  out  that  the  function  C has  been  derived  phenomenologically  for  the 

X *14'-16 

saturated  or  strong  scattering  case  by  Yura  and  Clifford  so  that  in 

principle  we  are  not  limited  to  the  first-order  expression  (62). 

We  generalize  (61)  to  the  nth  moment: 

<In(p)>  = n!  (o^r)  [ * • • [dp]  . . .dpTn  I"  II  <I(p\  ) >1  exp  ["  2 E C 1 

J J Li=l  i J L i^j  XijJ 


We  write  this  as 


(64) 


<l’ 

"(p)>  = n!<I(p)>nF 
n 

where 

= 1 

F2 

=Xfdp2dp4<I (p2) ><I  (^4)  > e4CX  ( 1 P2'P4 

1) 

[jdp2<I(p2)>]  2 

F 

i*  /•  __  ft 

J...J  dpjdp2,  . . dp"n  ^ n <I(p'1)>J 

exp 

[2  j,  VKV 

n 

[jdp<I(p)>]  n 

(65) 


(66) 


Our  assumption  on  C yields 


4Cx(|p2-P4l>  -1+4C (|p2-p4|) 

^ A 


(67) 


We  thus  simplify  F2  to 


4jJdp2dp4<I(p2)><I(p4)>Cv(|p2-plt  |) 
1+ — ^ 


[j  dp2<I(p2)>] 


(68) 


14. 

15. 

16. 


H.  T.  Yura,  "Physical  Model  for  Strong  Opticalr-Amplitude  Fluctuations 
in  a Turbulent  Medium",  JOSA,  64,  59,  January  1974, 

S.  F.  Clifford,  G.  R.  Ochs,  and  R,  S,  Lawrence,  "Saturation  of  Optical 
Scintillation  by  Strong  Turbulence",  JOSA,  64,  148,  February  3/  74, 

S.  F.  Clifford  and  H.  Tf  Yura,  "Equivalence  of  Two  Theories  of  strong 
Optical  Scintillation",  JOSA,  64,  1641,  December  1974. 
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For  the  denominator  for  the  focused  case,  we  write  from  Eq.  (10): 


:I(p>>  m(if  lu„l1 2  -§- J 


r2  / r_v»/3 

4a  2 ' ' 


(69) 


which  integrates  to 


1 


<I(p)>dp  = 2tt|U  I2  SSL 
o 2 


(70) 


The  integral  in  the  numerator  of  (68)  is  straightforward  but  laborious, 
We  state  here  the  result  for  F2 : 

/ Ian  ^5/3 

> 
o 


/us  \ 

L " 2a  2k2  - Vk 71 

F2  = l+(4ii)2k2  J uduj*  dse  0 $(u)sin2 

O 0 

where  the  approximation  relates  to  the  assumption  (67)  on  Cx . 
We  then  have  the  second  moment  from  Eq.  (65): 


(71) 


<I2(p)>  = 2<I(p)>2F2 
and  the  n ^ moment  is 

<i“>  = n! <I>n  jl+  (F2-l)  } 


(72) 


(73) 


Hence  nonunity  F2  represents  the  departure  from  an  exponential  intensity 
distribution. 

We  finally  derive  the  probability  distribution  using 

+ . . . 


< e-Sl>  - l-s<I>+  s2<l2>  s3<l3> 


2! 


3! 


(74) 


We  let  <I>  = Iq  and  find  the  probability  as 

“ T~ 


Pj(I)  - y-  e ° 
o 


1 + (F2-l)j"l-  + J-J  (75) 

L o o 

so  that  we  have  a first-order  correction  to  the  exponential  distribution. 
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We  note  that  F?  (and  F ) will  be  unity  for  weak  turbulence  and  for 

very  strong  turbulence  ( p ^-^0)  0 Also,  F2  will  be  unity  for  a small  focused 

source  (a  ),  which  physically  relates  to  a large  target  spot.  Conversely, 

0 

F2  will  depart  from  unity  and  therefore  show  the  effect  of  the  amplitude 
perturbation  term  for  the  case  of  intermediate  turbulence  and  a large 
source  (small  target  spot) . Physically  this  is  the  case  of  a quasi-point- 
source  attempting  to  scintillate  in  the  usual  manner  for  a point  source 
in  the  first-order  theory,5  but  nevertheless  interacting  with  the  speckles 
created  by  the  diffuse  target. 

The  exponential  distribution  applies  for  the  weak  turbulence  case 
simply  because  the  behavior  is  that  of  a diffuse  source  in  a vacuum 
(speckles).  It  applies  to  a large  target-spot,  implying  that  the  phase- 
perturbation-dominance  assumption  of  Section  II. A is  then  applicable,  i.e. 
that  the  target  (vacuum)  speckle  mechanism  dominates.  It  applies  for 
strong  turbulence  because  the  field  again  has  the  nature  of  that  from  a 
diffuse  source  (atmospheric  speckle  p^) . These  considerations  will  be 
clarified  below  and  in  Section  II. E. 

2.  Covariance  and  Variance  of  Intensity 

We  now  derive  the  general  covariance  function,  dropping  the  requirement 
that  C <<1. 

X /10N 

The  covariance  is  the  second  tern:  of  Eq.  (18): 

/v\4ff  _ _ rkpi-P2)  (P1-P2)  __ 

CjCPi.Pii)  “ \2ttl)  JJdp^dp4  <I(P2)>'I^P4^>  e H(p,p)(76) 


where 


P **  P2^P4  and  P = Pl"P2  • 

The  full  coherence  function  replacing  Eq.  (19)  is  (Appendix  A) 


H(p,p)  - 


p8/3_p8/3|  r 1 11 

P 5/3  Liwr  TFPT j 

o 


(77) 


2C  (p,p)+ 

A. 


2C  (p,-p) 

A. 
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where 


17 

1 

<yp,p)  = 0.132TT*k2Cn2L  j dt  j du  u~8/3sin2  [~~^1~t)] 

0 0 

•Jo  [ultP  + (l-t)p |]  (78) 

The  latter  expression  (78)  assumes  that  multiple  scattering  (saturation) 
does  not  apply  but  can  be  modified  for  such  a case,^  The  irradiance 
at  the  target  <I(p)>  is  given  by  (69)  for  the  focused  case  and  with  an 
additional  (kaQ/2L)2  term  in  the  exponent  for  the  collimated  case  as  before. 

We  again  omit  tedious  algebra  and  integrations  and  state  the  result 
for  the  focused  case: 


J*  r2dr2  H 

o 


ik 

— P-p 


H(p,p) 


(79) 


The  variance  is  given  by 


Cx(o)  = 


2r2^ 


p 5 /a 

o 


• | od“J„  ( t no)  e4Cx(p) 


These^general  results  involve  five-fold  integrals  (more  in  the  saturated 
case  ) and  have  not  been  numerically  evaluated.  A comparison  with  Eqs. 

(26)  and  (28)  show  that  the  amplitude  perturbation  term  simply  introduces 
the  terms  in  Eq.  (77),  i.e.  the  log  amplitude  covariance  for  a point 

17.  T. Wang, S.F. Clifford,  and  G.R.Ochs,  "Wind  and  Refractive  Turbulence  Sens- 
ing Using  Crossed  Laser  Beams",  Applied  Optics,  L3, 2602,  November  1974. 
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source.  A qualitative  interpretation  will  be  given  below. 

In  order  to  compare  these  results  with  those  of  Section  II.B.l,  we 
again  let  C^<<1,  and  we  find  that  (80)  and  (63)  yield 


o 


N 


1 + 0.528ir2k2C2L 
n 


du  u 


*'»  [f 


sin^ 


[ u2Lt(l-t) 

L 2k 


(81) 


This  is  consistent  with  Eq.  (72),  and  again  shows  the  departure  from  the 

unity  normalized  variance  obtained  when  C cannot  be  taken  as  zero. 

X 

3.  Further  Remarks  on  Probability  Distributions  and  Limiting  Conditions 


A systematic,  quantitative  description  of  parameter  realms  will  be 
given  in  Section  II. E.  However,  we  can  indicate  here  the  types  of  considera- 
tions involved.  Let  us  write  the  general  coherence  function  (77)  as 


H(p,p)  = e 0 H'(p,p) 


(82) 


where 


H 1 (p , p) 


(83) 


From  the  preceding  discussions,  we  know  that  non-negligible  terms  relate 
to  a nongaussian  field  (irradiance  is  not  exponentially  distributed,  and 
that  H,jsl  corresponds  to  a jointly  gaussian  field.  The  amplitude  perturba- 
tion effects  (C^)  will  be  influential  under  the  following  conditions  (see 
Eqs . (79)  and  (83)): 

h 

A.  Turbulence  is  moderate  or  weak  (i.e.  p >(L/k)  , and 

k 0 

B.  The  source  is  large  (aQ>(L/k)  )f 

These  conditions  can  be  reasoned  from  Eq.  (79)  by  considering  the  large 
range  of  rj  In  the  exponential,  leading  to  a small  range  of  p from  the 
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Bessell  function,  and  the  .orresponding  influence  of  C in  (83)  on  the 

x 

result  of  the  p integration.  Physically,  this  corresponds  to  nonsaturated 
turbulence,  with  a near-field  target  and  small  spot  acting  as  a quasi- 
point-source . 

Conversely,  if  is  small,  the  range  of  T2  is  limited,  and  that  of 

p expanded  such  that  does  not  make  an  important  contribution  to  the 

p integral.  This  means  that  a large  target  spot  results  in  dominance  of 

the  phase  perturbation  terms.  In  fact,  if  p »u  , H'  can  be  taken  as 

0 0 

unity,  the  field  is  jointly  gaussian,  and  the  covariance  is  given  by  Eq.  (47). 
Finally,  if  turbulence  is  strong  (Pq<  (L/k)*5)  , which  corresponds  to  the 
multiple  scattering  realm,  then  with  a large  or  small  source  the  o-integral 
will  be  controlled  by  the  phase  terms  in  (83)  . Thus  the  phase  perturbation 
dominates  in  strong  scattering;  the  field  will  be  jointly  gaussian  except 
in  the  transition  region  pq  £ oq,  where  it  is  simply  gaussian. 

4.  Further  Discussion  of  the  Covariance  Function  and  Scales 

The  qualitative  behavior  of  the  covariance  curve,  with  the  attendant 
covariance  scales  or  scintillation  patch  size,  can  be  described  as  follows. 
Suppose  that  we  have  a large,  focused  transmitter  (near-field  target), 
with  weak  turbulence  («Q .P^d/k)*5)  . The  scintillations  will  be  essentially 
those  of  a point  source  (unsaturated)  : 5 


0.124  C 2k//6  L11/6 


(84) 


for  the  log  amplitude  variance.  The  covariance  curve  for  irradiance  will 
be  as^shown  in  Figure  3a,  where  the  covariance  scale  is  proportional  to 
(L/k)  . Now  let  the  turbulence  strength  increase  (p  decrease)  until  p 
is  comparable  to  (L/k)*5;  the  log  amplitude  variance  peaks  at  approximately 
0.6  before  saturating,5  and  a new  covariance  scale  (Pq)  appears  as  concep- 
tualized in  Figure  3b.  Finally,  let  the  turbulence  become  stronger  yet; 
the  normalized  irradiance  variance  is  unity  and  the  covariance  curve  is 
as  given  in  Figure  3c.  The  third  possible  covariance  scale,  o itself 

1 G * 

Will  only  predominate  when  aQ  is  small  (far-field  target).  In  the  latter 
case  only,  the  speckles  will  not  be  transported  by  the  transverse  wind. 
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Figure  3.  Conceptual  normalized  function  of  irradiance  for  a large,  focused 
source  and  varying  turbulence  strength,  a 2 denotes  the  log 
amplitude  variance  for  a point  source.  * 

A.  Weak  turbulence  (p0>>(L/k)^) 

B.  Transition  region  (,p0%(L/k)^  . 

C.  Strong  turbulence  or  saturation  (p  <<(L/k)  . 


This  will  also  be  true  for  a large  target  spot  arising  from  defocusing  or 
collimat ion,  such  that  the  associated  speckles  are  small, 

II. C. Target  Glint? 

As  a preliminary  examination  of  the  effects  of  target  structure,  we 
present  in  this  section  an  analysis  of  the  mean  receiver  irradiance  in  the 
presence  of  ’’glints”  (specular  reflectors),  We  assume  that  the  target  is 
otherwise  diffuse,  and  that  phase  perturbations  predominate.  The  transmit- 
ter is  again  a TEM^  focused  or  collimated  source, 

We  again  start  with  the  field  at  the  receiver  as  described  by  Eqs.  (1,2). 
The  field  reflected  at  the  target  is  represented  as 


“<»>  - “diffuse  M + “glint  (P)-  Vp)+ 


where 


M 

U=  E U ^ (p  ) a exp 

glint  n incident  Mm  m 

m=l 


-(p-p)2 

m 


Ap 


m 


p - position  of  m^  giint 
m 

Ap  ” width  of  m glint 

® th 

a = complex  strength  of  m glint 
m 


Then  the  irradiance  is 


(85) 


W ■ ( 2Vx;)lf[Ud <0 i ) -HJg CP  1 )]  [u*(p2)+U?(p2)] 

• exp{g  [(Pi-p)2-(p2-p)2]  + i(i2  (p^p)^*  (p2,p)  ^ dpltdp2 

* (27L)2/iI“d^l)+“g(P.>]  [u3(P-2)+U.(p2,] 

. exp  1 [p2-p1  “2(pj-p2)pJ  + if2  (p2,p)+V2  (p2 .p) ^ dp: ,dp2 (86) 


We  then  simplify  the  field  at  the  target: 
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[ud(pi  )+U  (pi)]  [ud(P2)+Ug(P2)] 


* — 


* — 


<Ud(pl)Ud(p2)>  4 (Pl)U  (p2)>  + <U  (Pl)U  (p2)> 


o 8 


(87) 


Since 


* — 


<Ud(Pl)Ug(p2)>  = <ud(F1)><ug(p2)>  = o 

<Ug(Pl)ud(P2)>  = <Ug  (p"i ) ><U*  (p"2)  > = 0 
<ud(p’l)U*(p2)>  = <I(Pi)>  6 (pj-p2) 


^2(pi »p)+^*(F2,p) 

< e >= 


^ipl-?2.1)5/3 


(88) 


we  have 

<I(F)>  - <!(?)>„  + <i(F»g  (M) 

"here  <I(p)>d  Is  the  diffuse  term  analyzed  in  Section  II, A and  <I(F)>  ls 
the  glint  term:  8 


target : 


<I(p)>g~  ( 2ttl)  Wdpldp2  <Ug(Pi)U*(P2)> 

• exp  {g  [p^2(Pl-F2).p]-  ^ } 

p0  5/3  ’ (90) 

We  now  describe  the  complex  correlation  of  the  glint  field 


at  the 


* _ 


<Ug(0l)Ug (F2)>  - I la  |2  <Ttri,2„.rl  ‘>i+Pr-2g„-(pl'fp2H2P 

”'1  " ( A02 

m 


M 

+ £ a a <U(P  )U*(cT  )> 

mi?hn2  ml  m2  V m2; 
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► exp 


Pl~2p  *Pl+P  ^ P2~2p  . py+p  2 

ml  ®l  m2  z ra2 


Ap_ 


Ap_ 


(S 


m2 


where  for  the  focused  case 


-(if  l“ol2  t\  e V'U 


5/3 


(9 


The  corresponding  expression  for  the  collimated  case  has  an  additional 
(kaQ/2L)2  term  in  the  exponent.  The  coherence  function  can  be  readily 


shown  to  be 


<U(p  )U*(p  ): 

Tu  J m2 


> = k(r)2  iu0|2  ^ 2 L lpmi  Pra2^ ' 2L  ^Pmi  Pm* 


■ i 


-it  - Ik  r (pm  +p„  ) 


4a2  2L 


raj  m2 


dr  e 


F(r,p  -p  ) 
mj  m2 


(9: 


<U(VU*V>  ■ iv(r)2  I“„I2 


“o2  f - ' 4«2  "[2  L lr  (pmi 


dr  e 


£ r (v%)  - f (c4-%2) 


F(r,p  -p  ) 
m j m2 


(93 


where  the  mutual  coherence  function  F is 


2.91  C 2 L k2 


[ | Ip8^^8/3I  j 


F(r»P  -Pm  ) 
m j m2 


F(r,p)  = e 


(94! 


Finally,  the  mean  receiver  irradiance  (diffraction  pattern)  due  to  glints 
is  expressed  from  (90,91)  as 


M 


<xW>g 


V2iTlJ 


M 

E 

m2=l 


E 

m2=l 


a a*  <U(p 
nij  m2  mj 


)U*(p  )> 

IU2 


dPldp2  exp||^  [pi%^2(Pl-p2>p] 

% Apm2  \ p0  / 

where  clKp^ )U*(Pm2) > is  given  by  Eqs.  (93).  We  now  show  specific  cases 

1.  Single  glint  (M  = 1) 

From  Eq . (95)  we  have 


2 2 _ _ _ 
Pl'+p2_2p  .(pj+p2)+2p  2 

2. m 

Ap  2 
m 


With  the  following  change  of  variables: 


Pl-P2  *■=  p ) 2 2 I , 

> Pl+P2  = 2'(p2+4R2) 

P1+P2  = 2R  ) 


we  have 


The  first  integration  yields 


( 


(S 


! 


f RJo  Or?-  sr?lR] 

J m 


2tt 


-i&j 

v m ' 


Ap  2 
m 


dR 


Ap 


m 


-)•  ■ &■)  i 


Ap_ 


ZTT  -7  e 


m 


ik  — -- 
+ T—  Pm*  P 

L m 


(99) 


which  with  (98)  gives 


<w>  .(£)2 


<I(pm)>laml 


Ap  2 

m 


J 


f 


p j 


[r  lp-°mlp] 


1 r 

4 Ap_ 


Ap, 


m 


m 


(r)2]p2-(t) 


5/3 

dp 


(100) 


This  may  be  usefully  normalized  by  setting  p = p 


m 


The  result  (100)  may  be  interpreted  as  follows.  The  maximum  occurs 


at  p 


p^,  as  expected.  The  scale  of  the  diffraction  pattern  ( p - p |) 


is  seen  from  the  Bessel  term  to  be  L/kp  , where  p is  determined  from 

max  max 


the  exponential  as  the  smallest  of  Ap  , L/kAp  , and  p respectively.  The 

m m o 

corresponding  I(p)  scales  are  L/kAp  (diffraction  from  the  glint),  Ap 


m 


m 


(geometric  reflection  term),  and  p (effective  coherent  glint  size  in 

o 


strong  turbulence);  the  largest  will  predominate  in  each  case. 


2.  Arbitrary  Number  of  Glints  (M) 

This  involves  lengthy  but  straightforward  calculation.  The  result 


is 


M 


<I(P)> 


g 


l <I(p)> 
i=l 


S,+  * 


M 
E F 


i.J 


(I) 


(101) 


where 
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Fi,J(I) 


(r)2  R«[«'«'m1 


)U*(P  )>a  a (Ap  2 ) 

m m m r -M 


mj  mi  mj 


ij+' 


. exp<r- 


lpm  -p  I 

j ml  I 

Ap  2+Ap  2 
mi  "j 


I 


pdp  J (|y|p)  exp  | - j-  ■— 
° | 4 Ap 


“2  + Ap112+ 
ij+  1J  + 


and 


Ap 


Ap  2.A p 2 
m . 

— i i 

Ap  2+Ap  2 
ml~  m2 


(102) 


(103) 


Each  pair  of  glints  results  in  a cross-term  bright  spot  in  the  receiver 
irra^ance  locate  at 


This  is  illustrated  in  Figure  4. 
II. D. Incoherent  Illumination 


(104) 


We  now  depart  from  the  eeee  of  TE«0()  leeer  illumination  and  briefly 

multi  'I  V n0n"m°nOChr0matic  casa<  “hith  might  include  that  of  a highly 
multimode  Inner.  In  principle,  all  of  the  enpreanlonn  of  the  pre- 

ceding  nectionn  can  be  generalized  with  the  utilization  of  a 
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Irradiance 


Fourier  transform  or  partial-coherence  representation  of  the  source.  As 
the  source  becomes  monochromatic,  the  first-order  effect  will  be  to  wash 
out  the  speckles  (aQ  and  pQ  terms  in  Eq.  47).  This  will  obviate  the  influence 
of  the  phase-perturbation  terms  (Secs.  II. A, B)  and  the  remaining  scintil- 
lations will  be  due  to  the  amplitude-perturbation  term  and  to  geometrical 
considerations,  as  will  be  seen  below.  As  the  spectral  width  of  the  source 
grows  further,  scintillations  e.g.  for  a point  source  will  be  further 
smoothed  by  "spectral  diversity",  but  this  is  a weak  effect  requiring  a 
substantial  fractional  spectral  width. 

The  situation  may  be  simply  analyzed  if  we  assume  a completely  in 
coherent  but  quasi-monochromatic  target  spot.  We  start  with  the  intensity 
correlation  function  at  the  receiver  from  Eq.  (14),  and  utilize  uelta 
function  correlation  of  U(p})  with  U*(p2)  and  also  U(p3)  with  U*(p4)  to 
write  immediately 


<I(Pl)I(p2)>  - (^)jjdpidp3  |U(Pi')  |2|U(p3 


) 


. exp  | 4 Cx(p3~Pi >P)| 


(105) 


As  expected,  tne  target-speckle  (ctQ)  and  atmospheric  speckle  (pQ)  terms 
disappear,  leaving  the  log  amplitude  covariance  function.  Since  this 
function  is  unevaluated  within  the  Huygens-Fresnel  approach,  recourse  is 
again  necessary  to  first-order  or  saturation  theories  as  discussed  in 
Section  II. B.  It  is  interesting  to  note  however  that  Eq.  (105)  is  con- 
sistent with  Eq.  (8)  of  Ref.  19,  which  was  obtained  through  the  more 

. 13 

complicated  phase-screen  approach. 

In  the  case  of  a small  target  spot  (diameter< (L/k)  ) and  unsaturated 
scintillations,  Eq.  (105)  is  a statement  of  point-source  behavior.  For  a 
larger  target  spot,  it  appears  that  the  covariance  scale  will  be  smeared 
by  geometrical  effects,  and  this  can  be  directly  deduced  from  the  spatial- 
spectrum  representation  of  Ref.  19,  Eq.  (9).  (Note  however  that  the  quan- 

18.  D.L. Fried,  "Spectral  and  Angular  Covariance  of  Scintillation  for  Propa- 
gation in  a Randomly  Inhomogeneous  Medium" .Applied  Optics ,K), 721 , Apr . 1971 

19.  S.F.Clifford.G.R. Ochs, and  Ting-i  Wang, "Optical  Wind  Sensing  by  Observing 
the  Scintillations  of  a Random  Scene" , Applied  Optics ,14, 2844, Dec.  1975. 
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tity  q in  that  reference  is  necessarily  real,  contrary  to  the  notation 
there.  Also,  q and  y in  that  reference  may  be  identified  with  |u|  or 
I so  that  the  calculated  covariance  is  closely  related  to  but  not 
identical  to  that  for  I.)  It  is  known  from  other  analyses20  that  the 
scintillation  variance  is  reduced  or  "smoothed"  by  a featureless,  extended 
incoherent  source,  specifically 


' [“/tt/nj  '7/s  Hc 

where  D is  the  source  diameter.  This  is  also  the  "aperture  smoothing 
factor"  for  an  incoherent  receiver  in  unsaturated  scintillations.20 

In  the  case  of  saturated  scintillations  (po<(L/k)S,  the  covariance 
function  in  Eq.  105  must  be  suitably  modified.  The  behavior  of  the 
irradiance  correlation  function,  in  particular  its  scale  size,  is  not  well 

understood  for  this  case  of  incoherent  illumination  and  strong  turbulence 
scattering . 


II. E. Comparison  of  Alternative  Approaches; 
Behavior  vs.  Parameter  Realms 


Variance  and  Covariance 


The  basic  method  of  analysis  of  this  report  is  the  extended  Huygens- 
Fresnel  principle,  which  is  uniquely  suited  to  expressing  problems  of 
arbitrary  sources  in  terms  of  point-source  quantities.  It  also  has  the 
virtue  of  inherently  retaining  the  "saturation"  or  multiple  scat"  -ring 
effect  including  the  explicit  influence  of  the  quantity  p e.g  :or 
an  extended,  coherent  source. 

The  most  elegant  alternative  method  is  the  phase-screen  approach,13*19 
a first-order  (single-scattering)  analysis  which  also  handles  an  arbitrary 
source  (and  receiver)  geometry  but  with  severe  limitations  on  the  strength 
of  scattering  or  phas-  distortion  for  extended  coherent  sources.  It  appears 
that  this  approach  can  be  made  to  handle  saturation  effects  properly 
through  modification  of  a log  amplitude  covariance  function,14'17  although 

general  applicability  e.g.  to  p^effects  over  a coherent  source  has  not  yet 
been  clearly  demonstrated. 

20.  R.  F.  Lutomirski,  et  al,  Degradation  of  Laser  Systems  by  Atmospheric 
Turbulence,  Rand  Report  R-1171-ARPA/RC,  June  1973. 
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The  two  approaches  are  basically  complementary  when  properly  com- 
pared. The  Huygens-Fresnel  approach  is  much  simpler  in  its  handling  of 
extended  monochromatic  sources  (including  the  diffuse  target  case) , and 
also  simpler  in  its  setting  up  of  the  extended,  incoherent  case.  However, 
in  general  it  leaves  unevaluated  the  point— source  log— amplitude  covariance, 
which  must  then  be  obtained  from  the  first-order  theory  and  saturation 
arguments  described  above.  Also,  source  structure  may  be  advantageously 
represented  in  terms  of  its  spatial  spectrum,  as  shown  in  Ref.  19. 

A very  useful  "hybrid"  approach  would  be  to  utilize  the  saturation 

Cx  function  in  the  Huygens-Fresnel  results  for  various  sources.  Also, 

it  is  of  interest  to  explore  the  behavior  of  that  function  for  ncnmonochro- 

matic  sources:  how  sensitive  is  the  influence  of  p to  spectral  width, 

o 9 

i.e.  is  pQ  "washed  out"  as  readily  as  interferometric  phenomena  such  as 
speckles? 

We  are  now  in  a position  to  summarize  the  variance  and  covariance 
behavior,  for  a variety  of  sources,  in  a particularly  revealing  manner. 

We  will  identify  the  fundamental  independent  variables  and  all  possible 
permutations  of  their  mutual  inequalities  (asymptotic  parameter  realms) , 
and  attempt  to  show  the  limiting  value  of  the  variances  and  covariance 
scales  for  each  such  realm,  and  for  each  source  of  interest.  In  order  to 
do  this  we  will  utilize  results  from  phase-screen  or  first-order  theory 
and  from  Huygens-Fresnel  analyses  when  respectively  appropriate.  This 
type  of  summary  has  the  advantage  of  clear  physical  interpretation  and 
lacks  only  the  quantitative  details  of  the  transition  between  parameter 
realms.  Also,  in  some  cases  it  permits  the  deduction  of  e.g.  covariance 
scales  through  dimensional  analysis  or  the  requirement  of  consistency  with 
adjacent  parameter  realms. 

We  present  these  results  in  tabular  form,  with  associated  remarks 
below.  In  all  cases,  the  right  hand  column  represents  the  saturated  or 
multiple  scattering  realm  as  opposed  to  the  single-scattering  realm  in 
the  left-hand  column.  Also,  "near— field"  and  "far— field"  refer  to  the 
position  of  the  (point)  on-axis  receiver  relative  to  the  source  size.  The 
equalities  shown  omit  numerical  constants. 
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The  symbols  are  defined  as  follows: 

Path  length: 

Wavenumber : 

Coherence  radius: 

Gauss ian-beam  transmitter  radius: 

Incoherent  transmitter  diameter: 

Covariance  scale  (scintillation  patch 
size)  at  receiver; 

Beam  size  at  receiver: 


L 

k 


w 


Log  amplitude  variance: 

Normalized  irradiance  variance: 

Note:  for  log  normal  scintillations  oy 

Table  I:  XN 


Focused  TEMqq  beam.  Realm  (a)  is  the  first-order  point-source  realm, 

while  (b)  is  the  "transmitter  aperture  smoothing"  realm  within  the  first- 
24-27 

order  theory.  It  is  interesting  to  note  that  the  latter  has  the  same 

predicted  power-law  dependence  of  the  variance  on  the  transmitter  size  as 

that  for  the  incoherent  source  discussed  in  Section  II. D.  The  reason  for 

this  is  that  in  both  cases,  interference  does  not  play  a role:  in  the 

incoherent  case,  there  is  no  interference,  and  in  the  coherent  case,  with 

Pq>01q»  there  is  insignificant  wavefront  distortion.  Realm  (c)  represents 

significant  wavefront  interference,  a proliferation  of  diffraction  scale 
26 

patches,  and  beamspread  determiped  by  diffraction  from  p . Realms  (d,e) 
represent  a saturated  point  source,  ^ and  realm  (f)  represents  a large 
source  with  multiple  scattering.  In  the  latter  case,  which  has  not  been 
experimentally  verified,  the  covariance  scale  is  deduced  from  adjacent 
realms  (c,e)  by  the  requirement  of  consistency  at  their  respective  transi- 
tions. In  all  cases  involving  beam  breakup  or  saturation,  the  normalized 
variance  is  unity  in  agreement  with  a model  of  random-phase,  identical- 

g 

frequency  summed  oscillators. 

24.  D.L. Fried  and  J. B. Seidman, "Laser  Beam  Scintillation  in  the  Atmosphere", 
JOSA,  57^,  181,  February  1967. 

25.  A. Ishimaru,  "Fluctuations  of  a Beam  Wave  Propagating  Through  r Locally 
Homogeneous  Medium",  Radio  Science,  4^,  295,  April  1969. 

26.  J.R.Kerr  and  R.Eiss,  "Transmitter-Size  and  Focus  Effects  on  Scintilla- 
tions", JOSA,  62,  682,  May  1972. 

27.  J.R.Kerr  and  J.R.Dunphy,  "Experimental  Effects  of  Finite  Transmitter 
Apertures  on  Scintillations",  JOSA,  63_,  1,  January  1973. 
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Table  II: 


Collimated  TEM^q  beam.  Realms  (a,d,e)  are  identical  to  those  in  Table  I. 
Realms  (b,c)  are  essentially  plane-wave  illumination  within  the  first- 
order  theory,  while  (f)  represents  plane-wave  illumination  with  saturation. 

Table  III: 

This  is  the  case  analyzed  in  Sections  II. A, B,  in  which  a focused  TEMQ0 

beam  illuminates  a diffuse  target.  The  target  is  the  "source"  and  the 

receiver  is  in  the  transmitter  plane  (Figure  1).  "Near  and  far  field"  refer 

to  the  receiver-plane  relationship  to  the  illuminated  spot  on  the  target, 

which  in  some  realms  is  in  turn  determined  by  the  p diffraction  scale. 

o 

Realms  (a,b)  are  the  first-order  point  source  case  with  log  normal  scintil- 
lation statistics.  In  realms  (d,e)  the  "atmospheric  speckle"  predominates, 
with  joint  gaussian  field  statistics,  while  in  (c,f)  the  target  speckle 
predominates,  also  with  joint  gaussian  statistics.  In  realms  (c-f)  the 
model  of  independent  additive  oscillators  applies,  resulting  in  a normalized 

g 

variance  of  unity.  In  the  transition  region  between  (e,f),  the  field 
statistics  are  gaussian  but  not  joint  gaussian.  In  the  transition  region 
between  (a,d),  the  normalized  irradiance  variance  may  exceed  unity.  The 
amplitude  perturbation  term  (C  ) is  important  (Section  II. B)  in  realms 

A 

(a,b).  In  fact,  Figures  3a, c correspond  to  realms  a (or  b)  and  d,  respec- 
tively. 

Table  IV: 

Collimated  TEMqq  beam  illuminating  a diffuse  target.  This  case  is 
only  of  interest  for  oi^  >L/k,  such  that  "collimation"  has  meaning  (near 
field  conditions).  The  basic  considerations  are  the  same  as  for  Table  III. 

(We  also  note  here  a further  possibility,  that  of  a direct  monochromatic 
diffuse  source  (such  as  an  extended  laser  transmitter  with  a diffusing 
screen  in  front  of  it.)  In  this  case,  Table  I applies  with  the  only 

exception  that  w = «>  and,  in  realm  (b)  , o 2 = 1.) 

XN 

Table  V: 

Incoherent  (Nonmonochromatic)  Source.  In  this  case  we  assume  that 
the  spectral  width  of  the  source  is  not  sufficient  for  "spectral  diversity" 
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TABLE  V,  Incoherent  Source  (b) 
(Note:  w = °°) 


*/L/k<Pn  ^L7k>D 

o o 


8 


scintillations,  i.e.  smoothing  owing  to  broadband  effects.  Unfortunately, 

the  incoherent  case  is  not  yet  perfectly  understood,  Realm  (a)  is  a point 

source  in  the  nonsaturated  case.  Realms  (b#c)  correspond  to  an  extended 

source  also  in  the  single-scattering  case,  for  which  the  analysis  of  Ref, 

20  predicts  a -7/3  power-law  behavior  of  the  variance,  and  with  a covariance 

scale  that  is  related  to  the  source  size  through  geometrical  effects  (Ref, 

19  and  Section  II. D);  since  the  source  is  incoherent,  the  condition  p <b 

o 

in  going  from  realm  (b)  to  (c)  has  no  effect.  Realms  (d,e)  are  for  the 
saturated  point-source,  and  it  is  expected  that  as  the  spectral  width 
increases  over  that  for  the  monochromatic  case,  the  ’'atmospheric  speckles" 
(pQ)  are  washed  out,  although  the  degree  of  spectral  spread  required  to 
cause  this  is  not  understood.  The  variance  will  decrease  from  unity  when 
the  detector  has  insufficient  bandwidth  to  detect  the  "beats"  between  the 
individual  oscillators  comprising  the  nonmonochromatic  source.  Further 
collapse  is  expected  for  the  multiple-scattering,  near-field  case  of  realm 
(f).  As  mentioned  previously,  it  would  be  desirable  to  explore  (d-f)  fur- 
ther using  the  "saturation  covariance"  work  of  Refs.  14rl7. 

II. F. Further  Work  and  Applications 

The  application  of  the  efforts  described  in  this  report  is  to  the 
performance  of  coherent  optical  adaptive  (COAT)  systems.  The  desirable 
analytical  and  experimentally  verified  results  include  the  strength,  spatial 
scale,  and  spatial  and  temporal  spectrum  of  scintillations  and  their  sta- 
tistics, for  a variety  of  sources  including  single-  and  multi-mode  illumi- 
nated diffuse  and  structured  targets.  The  coherent  or  complex  amplitude 
correlation  (mutual  coherence  function)  of  the  receiver  field  is  equally 
pertinent . 

Finally,  the  turbulence  effects  should  be  applied  to  control  system 
models  of  actual  adaptive  systems  which  are  generally  phase  correcting  sys- 
tems. We  have  given  this  preliminary  attention  with  the  assistance  of  sys- 
28 

terns  engineers,  including  the  dynamics  of  target  lock-on  and  the  resultant 
evolution  in  spot  and  speckle  conditions, 

Significant  progress  has  been  made  on  each  of  these  problems,  as 
described  above, 

28.  We  acknowledge  witn  appreciation  the  collaboration  of  James  E,  Pearson 
and  others  at  the  Hughes  Research  Laboratories,  Malibu^  California, 
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We  detail  below  specific  topics  requiring  further  analytical  work: 
Completion  of  Fundamental  Understandings 

A.  Complete  the  coherent-source,  diffuse  target  covariance  analysis 
of  Section  II. B.,  utilizing  the  full  point-source  mutual  coherence  function 
including  terms.  This  primarily  reduces  to  thorough  numerical  analysis 
of  the  integral  results;  also,  the  "saturation"  form  of  C (Refs.  14-17) 
should  be  incorporated.  The  results  should  be  extended  t o the  time-lapsed 

covariance,  temporal  and  spatial  spectra,  and  probability  distribution  of 
irradiance. 

B.  Complete  a similar  analysis  with  glints  and  other  target  structure, 
i.e.  through  second-irradiance  moments  and  spectra  (Section  II.C.),  and 
including  the  complex  amplitude  (mutual  coherence  function) . 

C.  Incorporate  incoherent  receiver-aperture  smoothing. 

D.  Incorporate  multimode  sources,  and  the  extreme  of  incoherent  sources 
(Section  II. D.)  including  the  saturation  form  of  C . In  particular,  the 
incoherence  or  spectral  width  required  to  wash  out  "atmospheric  speckles" 

(Pq  scintillation  scale)  should  be  understood. 

Applications 

A.  Characterize  target  models  at  wavelengths  of  interest  and  incorporate 
them  into  the  analysis. 

B.  Analyze  the  effects  of  moving  targets. 

c.  Analyze  the  influence  of  all  turbulence  effects  on  the  operation 
of  real  adaptive  systems,  including  lock-on  dynamics. 

Secondary  Aspects 

A.  Influence  of  <x4»  or  the  phase-amplitude  correlation  in  the  point- 
source  mutual  coherence  function, 

B.  Effects  of  inner  scale  of  turbulence. 

III.  Experimental  Effort 

The  experimental  effort  is  proceeding  in  three  stages.  The  first 
involves  preliminary  experiments  over  a short  (100m)  path  with  a Helium- 
Neon  laser  operating  at  6328  X,  with  a quasi-diffuse  (Scotchlight)  ta~,et. 
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This  project  was  designed  to  permit  shakedown  of  the  experimental  tech- 
niques and  system  electronics  and  verification  of  basic  theoretical  pre- 
dictions at  an  early  date,  and  has  been  completed.  The  second  ecage 
involves  operation  over  a 1 km  uniform  path  with  an  Argon  laser  operating 
at  4880&  using  a similar  target.  The  required  installations  are  nearly 
complete  at  our  instrumented  field  site.  The  third  stage  involves  the  use 
of  a TEA  laser  (3.5p  and  10. 6p)  which  is  currently  being  fabricated  and 
which  will  also  be  utilized  at  the  field  installation.  Characterization 
of  target  models  will  also  be  conducted  in  the  laboratory  at  wavelengths 
of  interest,  and  target  structure  will  be  added  at  the  field  si^e 

A basic  problem  in  these  experiments  is  the  energy  available  off  of  a 
diffuse  or  q aasi-dif fuse  target  at  a distance.  At  visible  wavelengths, 
the  directional  properties  of  Scotchlight  permit  an  adequate  signal-to- 
noise  ratio  while  preserving  the  pertinent  properties  of  a "diffuse” 
target;  that  is,  the  speckles  in  the  absence  of  turbulence  have  identical 
properties  to  those  observed  for  a truly  diffuse  source.  At  the  inrrared 
wavelengths,  the  high  instantaneous  peak  power  of  the  TEA  laser  is  required. 

The  project  goal  is  the  observation  and  measurement  of  a number  of 
statistical  properties  as  described  in  the  analytical  discussion  of  Section 
II.  The  basic  system  is  illustrated  in  Figure  5.  The  coherent  laser 
illumination  is  modulated,  focused,  spatially  filtered,  and  expanded  to 
impinge  on  the  remote  target  ~urface . The  output  beam  may  be  focused, 
collimated,  or  operated  in  an  intermediate  condition.  The  return  energy 
is  detected  by  two  separable  receivers  located  adjacent  to  the  transmitter, 
demodulated,  filtered,  and  recorded  on  analog  tape  for  digital  processing. 
Simultaneously,  the  correlation  of  the  two  irradiance  signals  is  computed 
electronically  with  variable  time  averaging  and  displayed.  The  measure- 
ments are  repeated  at  various  receiver  separations  to  yield  a covariance 
curve,  and  the  turbulence  strength  is  measured  using  a microthermal  prooe 
system.  The  output  from  final  processing  includes 

Turbulence  level  (C  2 and  p ) 

n o 

Spatial  covariance  function 

Probability  distribution  of  irradiance 


Modulator 


Normalized  variance  of  irradiance 
Temporal  power  spectrum 
Meteorological  parameters 

As  discussed  in  Section  II,  there  are  in  general  six  possible  parameter 
realms  resulting  from  variations  in  transmitter  configuration,  wavelength, 
range,  and  turbulence  strength.  It  is  our  intent  to  cover  all  of  these 

realms;  as  many  as  possible  will  be  covered  at  each  of  the  operating  wave- 
lengths . 

Initial  Results 

The  purpose  for  the  short-path  experiments  was  to  aid  in  the  design 
and  shakedown  of  the  optical,  electronic,  and  data  processing  systems  to 
be  used  in  the  further  stages.  The  low  path-integrated  turbulence  levels 
over  such  a path  in  the  winter  months  result  in  large  values  of  p compared 
to  transmitter  size  ( uq ) and  Fresnel  scale  (L/k)'4.  The  relevant  parameters 
are  listed  in  Table  VI,  and  it  is  noted  that  the  basic  parameter  realm  was 
limited  to  (L/k)  ^ aQ  < pq.  However,  defocusing  of  the  beam  permitted  the 
measurement  of  smaller  speckles. 

Visual  observations  indicate  the  expected  speckle  pattern,  and  show 
that  the  pattern  is  nearly  sationary.  This  indicates  that  the  "target 
speckles"  (ao)  are  predominating.  In  order  to  achieve  reasonable  time 
averaging,  it  therefore  proved  necessary  to  rotate  the  target  so  that  the 
speckle  field  is  swept  across  the  detectors.  The  rate  of  motion  is  limited 
by  the  receiver  bandwidth. 

Covariance  functions  of  irradiance  have  been  obtained  under  a number 
of  conditions.  Figure  6a  sho-s  the  measured  covariance  function  for  a 
focused  beam  over  the  85m  path  length,  with  a = 0.013m,  C 2 = 5.0xl0~15m2/3 
and  po  = 0.15m.  The  theoretical  curve  is  also  shown,  based  on  the  reasonings 
of  Section  II. A.  Since  the  point-source  log  amplitude  variance  (o  2)  was 
quite  low,  with  a comparable  size  for  the  target  speckles  and  Fresnel  scale, 
the  amplitude  perturbation  term  of  Section  II. B.  was  not  pertinent  (see 
Figure  3).  Figure  6b  shows  a similar  curve  with  a defocused  beam  such  that 
the  spot  radius  on  the  target  is  5mm  compared  to  1.3mm  in  the  preceding 
case.  The  expected  decrease  in  speckle  size  is  observed. 
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TABLE  VI.  Experimental  Parameters 


Wavelength 

Range 

Transmitter  Beam  Radius  a 
/L7k 

/XL 

Range  of  values 

Laser  Power 
Modulation  Frequency 
Signal  Bandwidth 
Receiver  Size 
SNR  (Scotchlight) 

SNR  (True  Diffuse) 


Stage  1 

Stage  2 

,6328p 

0.488p 

85m 

100- 1000m 

.013m 

.005-. 02m 

•0029m 

.0028- 

•0088m 

.0073m 

.007-. 022m 

0.K1.0 

meter 
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Figure  6a.  Measured  covariance  of  irradiance  for  focused  beam  and 
parameters  shown. 
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Figure  6b.  Measured  covariance  of  irradiance  for  defocused  beam  and 
parameters  shown , 
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Target  Plane 
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Figure  7.  Variables  for  mutual  coherence  function. 


Long- time-averaged  6 minutes)  measurements  of  the  normalized  variable 

of  irradiance  are  summarized  below,  for  the  focused  case* 


C 2(nf2/3) 
n 

a2/ 1: 

-1ft 

8.2x10 

1,12 

-1  ft 

8.9x10 

0,85 

4.5x10  ^ 

1.16 

The  theoretical  value  is  of  course  unity. 
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Appendix  A.  Generalized  Mutual  Coherence  Functions 

The  point-source  mutual  coherence  functions  are  crucial  to  the  utiliza- 
tion of  the  extended  Huygens-Fresnel  approach.  In  this  Appendix  we  review 
the  generalized,  high-order  mutual  coherence  function.  The  derivation  will 
not  be  given  but  is  based  upon  considerations  suggested  by  Yura. 

The  general  nth  order  mutual  coherence  function  is 

H(pi , P2 • • • P2n;Pl *P2 • • • P2n) 

= <exp  j*<PitPl)+  <I>*<P2.P2>+---+  *<P2n..1.P2n_i>+  **<p2n»P2n^  > (Al) 


where  = X + • 

We  now  make  the  following  assumptions: 

1.  x bas  a j°int  normal  distribution,  with  <x>  = "°x  * 

2.  <J>  has  a joint  normal  distribution,  with  «t>>  = 0. 

3.  x and  <f>  are  independent. 

There  are  limitations  on  the  validity  of  these  assumptions.  It  is  known 
that  x is  not  gaussian  in  the  case  of  saturated  scintillations,  and  as 
pointed  out  in  the  preceding  report  on  this  program,  it  is  in  general  not 
jointly  normal.3  Also,  <X*>  is  in  general  nonzero  and  can  be  calculated 

from  the  first-order  theory. 

Assumption  #1  can  be  replaced  with  the  assumption  of  a gaussian  refrac- 
tive index  field,  bu‘.  this  is  a very  poor  representation  of  the  actual 
nature  of  that  field.22  Alternatively,2  we  may  simply  say  that  the  results 
to  be  given  here  are  correct  to  "second  order  in  the  refractive  index  varia- 
tions", which  is  at  best  a heuristic  approximation  in  conditions  of  saturated 
scintillations.  However,  we  will  not  dwell  further  here  on  the  degree  of 
approximation  involved  in  the  assumptions,  which  comprises  a fruitful  area 
for  further  work.  It  is  sufficient  to  point  cut  that  very  useful  results 

are  obtained. 

The  primary  result  is  as  follows: 


21. 

22. 


H.  T.  Yura,  private 
J.R.  Kerr,  et  al,  " 
through  Atmospheric 


communication. 

Propagation  of  Multiwavelength  Laser  Radiation 
Turbulence",  RADC-TR-74-320 , November  1974, 


(A003340) . 
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'ij  j=l  i=l  X2j+i,i 


(A2) 


where  M is  the  number  of  j satisfying  2j+i=2n  (for  i even)  and  2j+i=2n-l 

(for  i odd).  The  wave  structure  function  and  log  amplitude  covariance 

function  C comprise  a generalization  of  those  normally  seen  in  the 
X23 

literature.  Referring  to  Figure  7,  they  refer  to  complex  phase  or  log 
amplitude  quantities  as  a function  of  (p^-p^)  the  observation  plane, 
for  fields  emanating  from  and  p^.  Specifically, 
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(A4) 


It  may  be  noted  that  the  functions  satisfy  (p,p)  reciprocity. 
Special  Cases 

(a)._  n_=i  -y 

H(Pi  ,P2;P:  .p,;)  = e ^12 


(A5) 


(b)  n = 2 

H(Pi ,p2,p3,p4;p1 ,p2,p3,p4) 


e“^s(Di  2~Dj  3+Dj  4+D23-D24+D34)+2C^  ^+2C 


23.  D.  L.  Fried,  "Atmospheric  Modulation  Noise  in  an  Optical 
Receiver", IEEE  Trans,  on  Quantum  Electronics,  QE-3,  213, 
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This  is  a generalization  (double-arguments)  of  the  result  of  Ref.  23 


H = e4CX  (P3~Pl) 


☆U.8.  GOVERNMENT  PRINTING  OFFICE1  1976-614-357/196 


(c)  n = 2:  Pi=P2»P3=Pit (calculation  of  <I1I2>) 

See  Eqs . (18,20)  in  Section  II. 

(d)  n = 2;  Pi=P2»P3=P4;Pi=P2=P3I=P4(calculation  of  f°r  point 

source) 
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